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We show that an atomic system in a periodically modulated optical trap displays an ideal mean- 
field symmetry-breaking transition. The symmetry is broken with respect to time translation by the 
modulation period. The transition is due to the interplay of the long-range interatomic interaction 
and nonequilibrium fluctuations. The observed critical dynamics, including anomalous fluctuations 
in the symmetry broken phase, are fully described by the proposed microscopic theory. 
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The mean-field approach has been instrumental for de- 
veloping an insiglit into symmetry breaking transitions in 
thermal equilibrium systems [1] . It has been broadly ap- 
plied also to nonequilibrium systems, the studies of pat- 
tern formation being an example [2]. However, close to 
the phase transition point the mean field approximation 
usually breaks down. This happens even in finite-size 
systems provided they are sufficiently large. 

In this paper we study a nonequilibrium system of 
~ lO'' particles which, as we show, displays an ideal 
mean-field symmetry breaking transition. It is accompa- 
nied by anomalous fluctuations, which are also described 
by a mean-field theory. The system is formed by moder- 
ately cold atoms in a magneto-optical trap (MOT) [3, 4]. 
The atoms are periodically modulated in time [5]. Peri- 
odically modulated systems form one of the most impor- 
tant classes of nonequilibrium systems, both conceptually 
and in terms of applications. They have discrete time- 
translation symmetry: they are invariant with respect to 
time translation by modulation period Tp. Nevertheless, 
they may have stable vibrational states with periods that 
are multiples of Tp, in particular 2Tp. Period doubling is 
well-known from parametric resonance, where a system 
has two identical vibrational states shifted in phase by 
TT. It is broadly used in classical and quantum optics and 
has attracted significant attention recently in the context 
of nano- and micro-mechanical resonators [6-9]. 

In a many-body system, dynamical period doubling in 
itself does not break the time translation symmetry. This 
is a consequence of fluctuations. Even though each vi- 
brational state has a lower symmetry, fluctuations make 
the states equally populated and the system as a whole 
remains symmetric. However, if as a result of the inter- 
action the state populations become different, the sym- 
metry is broken. It is reminiscent of the Ising transition 
where the interaction leads to preferred occupation of 
one of the two equivalent spin orientations, except that 
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the symmetry is broken in time. For atoms in a para- 
metrically modulated MOT spontaneous breaking of the 
discrete time-translation symmetry was observed exper- 
imentally by Kim et al. [10]. 

Here we show that the time-translation symmetry 

breaking in a modulated MOT results from the coop- 
eration and competition between the inter-particle inter- 
action and the fluctuations that lead to atom switching 
between the vibrational states. The interaction is weak. 
On its own, it cannot change the state populations. How- 
ever, it may become strong enough to change the rates 
of fluctuation induced switching, which in turn will cause 
the population change. We provide a quantitative theory 
of the phenomenon. Wc measure the critical exponents 
and the frequency dispersion of the susceptibility. The 
observations are consistent with the mean-fleld behavior 
and are in good agreement with the theory. 

In the experiment, ^^Rb atoms in a MOT were cooled 
down to «0.4 mK. Full three-dimensional conflnement 

was achieved with three pairs of counter-propagating 
laser beams. The intensity of the beams along the MOT 
axis z was 0.19 mW/cm^, the transverse beam intensities 
were 5 to 10 times larger. The magnetic field gradient at 
the trap center was 10 G/cm. The transverse beams were 
detuned from the atomic transition by (5 ~ — 2.3rp, the 
longitudinal beams were further detuned by 5 MHz (the 
atomic decay rate is rp/27r « 5.9 MHz). The atomic 
cloud motion was essentially one-dimensional, along pr- 
axis, with frequency a;o/27r w 45 Hz and the damping 
rate F « 36 s~^. The total number of trapped atoms A^tot 
was varied by changing the hyperfine-repumping laser in- 
tensity at a decrease rate of 0.5% per second. For the or- 
der parameter and the variance measurements very close 
to criticality, the decrease rate was reduced to 0.03%. 

The beam intensities were modulated at frequency 
cop = 2uJo by acousto-optic modulators. When the modu- 
lation amplitude exceeded a threshold value, after a tran- 
sient, the atomic cloud split into two clouds [5], which 
were vibrating in counter-phase with frccjuency ojp/2, see 
Fig. la. We took snapshots at the maximum separation 
of the clouds at the frame rate of 1-50 Hz to obtain the 
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population of each cloud. For a small total number of 
trapped atoms TVtot the populations were equal and the 
cloud vibrations had equal amplitudes. 
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FIG. 1: (Color) (a) The two-dimensional density profiles of 
the vibrating clouds for the maximal displacement along the 
MOT axis (« 4.5 mm). The lower, middle, and upper panels 
refer to the symmetric, close to critical, and broken-symmetry 
states, with iVtot « 1.5 x 10^, 5.2 x lO'', and 6.7 x 10® , respec- 
tively; Nc ~ 5.6 X 10®. The cloud centers are indicated by 
the ticks at the bottom of the panels and are sketched in the 
right panel, (b) The amplitude changes of the cloud vibra- 
tions as functions of the population of the other cloud. The 
red dashed lines show the expected linear dependence. The 
black dashed line shows A''i = = Nc/2. The amplitude 
of the cloud with smaller population beyond the transition, 
cloud 1, monotonically increases with A'tot, whereas that of 
the larger cloud, cloud 2, is maximal at criticality. The error 
bars show the standard deviation for 20 measurements. 

Once A^tot exceeded a critical value Nc {Nc ^ 10^ in 
our experiment), the average populations of the clouds 
became different, as observed earlier [10]. The vibra- 
tion amplitudes became different, too, see Fig. lb. This 
is spontaneous breaking of the discrete time-translation 
symmetry, the system becomes invariant with respect to 
time translation by 2tp, not the modulation period Tp. 

As phase transitions in equilibrium systems, the ob- 
served transition is a many-body effect. This is evidenced 
by the required critical number of atoms, which is finite 
but large, typical of phase transitions in cold atom sys- 
tems [11-15]. However, in contrast to systems in ther- 
mal equilibrium, our system is not characterized by the 
free energy and has a discrete time-translation symme- 
try. Most importantly, there are no long-range corre- 
lations of spatial density fluctuations. The intra-cloud 
density fluctuations are irrelevant, they are uncorrelated 
with the inter-cloud fluctuations, which, as we show, are 
responsible for the transition and display critical slowing 



down. As a result, the symmetry breaking is quantita- 
tively described by the mean-field theory. 

The irrelevance of intra-cloud density fluctuations is a 
consequence of the time scale separation and the weak- 
ness of the interaction. The decay of correlations of intr- 
acloud density fluctuations is characterized by the damp- 
ing rate F of atomic motion, which is determined by the 
Doppler effect in the MOT [4] and does not change at 
the phase transition. On the other hand, switching be- 
tween periodic states in modulated systems requires a 
large fluctuation and involves overcoming an effective ac- 
tivation barrier [8, 16-20]. For low temperature, switch- 
ing is a rare event. We measured the rates Wnm of inter- 
cloud n ^ m switching (n, m = 1, 2) directly for low iVtot 
from the decay time of the difference of the cloud pop- 
ulations, which gave W12 = W21 ~ 1 s~^. Therefore, in 
the experiment there is a strong inequality between the 
characteristic reciprocal times, uip ^ Wnm- 

During switching an atom is far from the clouds. Its 
interaction energy with individual atoms in the clouds 
is much smaller than kgT. However, the total energy of 
interaction with all atoms in the clouds may exceed fc^T, 
and then it changes the inter-cloud switching rate. The 
interaction depends on the density distribution within 
the clouds and ultimately on the cloud populations iVi 2. 

The switching rates modification by the atom-atom in- 
teraction is the underlying mechanism of the symmetry 
breaking. Since these are the cloud populations A^i,2 that 
control the switching rates and that are changed as a re- 
sult of switching, we choose their mean relative difference 
77 = {N2 — Ni)/Ntot as the order parameter. A natural 
control parameter is the reduced total number of atoms 
9 = (TVtot - Nc) /Nc. 

Our main observations are presented in Figs. 2 and 3. 
They show that all measured quantities display mean- 
field critical behavior. The order parameter scales as 
\ri\ cx 6"^ for > 0, with /? sa 1/2, see Fig. 2a, and 
?7 = for 6* < (this scaling was also seen earlier [10], 
but now we have come much closer to the critical point, 
and the precision is much higher). Fluctuations of the 
population difference are characterized by the variance 
CT^ = N^^l{{N2 - Nif) - if. It scales as cr cx |6l|-'' with 
7 « 1 for positive and negative 9, see Fig. 2b. Equi- 
librium systems at criticality show a strongly nonlinear 
response to the symmetry-breaking static field, like a 
magnetic field at the Ising transition. For our system 
an analogue of such field is an additional periodic driv- 
ing at half the strong-field frequency. With this field, 
the overall time-translation period is 2tp , and the vibra- 
tional states of period 2tp become non-equivalent [21]. 
We studied the effect of the dynamic symmetry breaking 
field at criticality, = 0, by asymmetrically modulat- 
ing the counter propagating beams. Fig. 2c shows that 
the cloud population difference scales with the additional 
field amplitude h as |?7| cx h^/^ with (5 ~ 3, in agreement 
with the mean field theory. 

The critical slowing down should lead to a strong re- 
sponse of the population difference to a field detuned by 



3 




FIG. 2: (Color) (a) The relative difference of the cloud pop- 
ulations ri — {N2 — Ni}/Ntot as a function of the control 
parameter 9 — (A^tot — Nc) /N^; the critical total number of 
atoms Nc — 1.6 x 10^. The solid curves show the mean- 
field theory, r] = tanh [(e -|- 1)77]. For < < 1 |r;| oc 6^ 
with P — 0.51 ± 0.01 (inset), (b) The variance of the order 
parameter = lO^cr'^ and its scaling (inset). The critical ex- 
ponents in the symmetric and broken-symmetry phases are, 
respectively, 7 = 1.04 ± 0.21 and 7 = 1.11 ± 0.13. (c) The 
order parameter at criticality {6 = 0) as a function of ampli- 
tude h of the additional modulation at frequency 01^/2; h is 
scaled by the strong modulation amplitude. The solid line is 
77I oc h^^^ , with S — 3.0 ± 0.8. (d) Order parameter \t]\ as a 
function of A'tot (in units of 10^ atoms) and the intensity of 
the additional resonant radiation (scaled by the saturation in- 
tensity, 3.78 mW/cm^). The red curve is the phase transition 
line. The error bars in panels (a) and (c) show the standard 
deviations of 10 measurements. 



a small frequency Q from ujf/2. This is an analog of ap- 
plying a weak slowly varying field to thermal equilibrium 
systems. The amplitude and phase of the forced oscilla- 
tions of r] at frequency fl as functions of 6 are shown in 
Fig. 3. They display characteristic asymmetric resonant 
features, again in excellent agreement with the mean- field 
theory. 

In general, one should not expect that our oscillating 
system can be described by the Landau-type theory with 
standard critical exponents, and in fact the theory has 
to be extended. We explain the observations as due to 
kinetic many-body effects. They turn out to be strong as 
a consequence of the exponential sensitivity of the inter- 
cloud switching rates Wnm- 

If the atom-atom coupling can be disregarded, the 
switching rates are equal by symmetry, = = 

= Cexp {-R'^^^/kBT) with C - Wd, where uj^ ~ 
3r is the frequency of small-amplitude damped atom 
vibrations about the cloud center. Of primary inter- 
est is the activation energy For weekly non- 
sinusoidal period-2 vibrations it was calculated and mea- 
sured for single-oscillator systems with cubic nonlinear- 
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FIG. 3: (Color online) The amplitude (a) and phase (b) of os- 
cillations of the order parameter induced by an extra modu- 
lation at frequency up/'i + Q with 0=0.1 Hz. The solid curves 
show the theory with the A'tot — >■ switching rate used as a 
fitting parameter. The error bars show the standard devia- 
tions of 100 measurements. 



ity [8, 16, 17]. If we use the same model, for the 
present experimental parameters the theory [16] gives 
i?(°V^B -3.4 mK. This is within a factor of 2 of the 
value extracted from the measured W'^'^^ and estimated 
T, which is satisfactory given that in the present case the 
parametric modulation was comparatively strong (the 
squared MOT eigenfrequency ujq was modulated with rel- 
ative amplitude 0.9), and therefore the vibrations were 
noticeably non-sinusoidal. 

The atom-atom interaction changes the switching ac- 
tivation energy. The change becomes significant only 
because of the cumulative effect of the interaction of a 
switching atom with the atoms in the clouds. For not 
too strong interaction the rate of 1 ^ 2 switching be- 
comes 



Wi2iNi; Ntot) = exp[a(7V2 - iVi) + PNu 



(1) 



with a,/? cx 1/T (to obtain W12 one should replace 
iVi o The parameters a,/3 are determined by 

the generalized work of the interatomic force done on 
an atom during its intercloud switching, see Supplemen- 
tary Material. The probability Pi to have iVi atoms 
in cloud 1 for given A'tot can be found from the bal- 
ance equation. For A'tot ^ 1, the stationary probabil- 
ity Pf*(A^i) is a function of the reduced population dif- 
ference X — {N2 — A^i) /A'tot and the control parameter 
9 = a A'tot — 1. It has a sharp maximum whose position 
gives the order parameter 77 = (x) . 

For fixed A'tot, near criticality (|x|, 16*1 <C 1) the station- 
ary distribution has a characteristic Landau mean-field 
form Pf cx exp [- A'tot (x" - 6ex^)/l2] [22]. For 6 < 
it has a peak at x — rj — 0, which corresponds to equal 

mean cloud populations. For 9 > and |6'| 3> A', 
the symmetry is broken and (x) ^ rj ~ ±{39)^^^ 
describes the dependence of on 9 in Fig. 2a. The form 
of Pf* explains also the scaling oc |6'|~^ of the variance 
cr2 = (a;2j _ 2b. 

The critical total number of atoms (where 6* = 0) is 
Nc = a^^. Thus, Nc oc T. This dependence was tested 
by heating up the atoms with additional resonant radi- 
ation. Indeed, Nc was found to increase almost linearly 
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with the radiation intensity, see Fig. 2d. 

An extra modulation of the beam intensities 
h cos{ujFt/2+(/)h) leads to factors exp(/ii2) and exp(— /ii2) 
in the switching rates W12 and W21, with hi2 oc h/T [21]. 
In calculating these factors one can disregard the weak 
atom- atom interaction. As a result, Pf* is multiplied by 
exp(Aftot/^i2a;). This gives r] oc h^^^ at criticality {9 = 0), 
in agreement with the experiment. If the frequency of 
the extra modulation is slightly detimcd from uopjl, the 
linear response does not diverge at criticality. The result 
of the mean-field calculation in the Supplementary Ma- 
terial is in excellent agreement with the experiment, as 
seen from Fig. 3. 

An independent estimate of the interaction strength 
can be obtained from the dependence of the vibration 
amplitudes of the clouds on the cloud populations. For 
comparatively weak interaction, the interaction-induced 
change of the amplitude of an nth cloud should be pro- 
portional to the number of atoms in the other cloud N^-n 
{n = 1,2). This is indeed seen in Fig. Ic. The main 
contribution comes from the long-range interaction, the 
shadow effect [4, 23, 24]. This interaction also gives the 
main contribution to the parameters a, /3 and thus deter- 
mines the critical total number of atoms Nc- The value of 
Nc obtained in a simple one-dimensional model [16, 21] 
that assumes sinusoidal vibrations is within a factor of 2 
from the experimental data, as are also the slopes of the 
straight lines in Fig. Ic. 

In the experiment, the total number of trapped atoms 
was slowly fluctuating. The mean-field theory still re- 
mains applicable, but needs to be extended, sec Sup- 
plementary Material. An important consequence is that 

— 1/2 

in the low-symmetry phase, 1 ^ 6 ^ Nc , the vari- 



ance is modified, NcO-"^ = (29)-'^ + 3(49 + e)~'^ , where 
e = Wout/W^^°' exp(/37Vc), e < 1, with Wont being the 
probability for an atom to leave the trap per unit time. 
The variance is larger than in the symmetric phase for 
the same 16*1 by factor 5/4 for ^ » e and is smaller by 
factor 2 for smaller 9, i.e., closer to Nc- The scaling 
cr^ oc 9~^ holds in the both limits. This is the scaling 
seen in Fig. 2b; we could not come close enough to A^c to 
observe the crossover; in the broken-symmetry phase 
remained larger than in the symmetric phase. 



In conclusion, we demonstrate the occurrence of an 
ideal mean-field transition far from equilibrium, the 
spontaneous breaking of the discrete time-translation 
symmetry in a system of periodically modulated trapped 
atoms. The mean-field behavior is evidenced by the crit- 
ical exponents of the order parameter and its variance, 
as well as the nonlinear resonant response at criticality 
and the susceptibility as a function of the distance to the 
critical point. We explain the effect as resulting from the 
interplay of the interaction and nonequilibrium fluctua- 
tions, where the interaction, even though comparatively 
weak, is strong enough to affect the rates of fluctuation- 
induced transitions between coexisting vibrational states. 
The proposed theory is in full agreement with the obser- 
vations. 
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PACS numbers: 

I. THEORY OF INTER-CLOUD SWITCHING 

A. Single-atom switching 

Here we provide a detailed closed-form description of 
the critical phenomena in a parametrically modulated 
MOT. We start with the single-atom dynamics for reso- 
nant MOT modulation, where the modulation frequency 
Wi? is close to twice the eigenfrequency Wq of atomic vi- 
brations in the MOT. The modulation can lead to the 
onset of period- two vibrations at frequency Wir/2. These 
vibrations are close to sinusoidal for moderately strong 
modulation, there is no dynamical chaos in the range of 
phase space of interest. In the single- atom approximation 
the equation of motion can be written as 

mat = G(r,r;i) + f(f) (1) 

where r is the position vector of an atom, rua is the 
atomic mass, and G is the well-understood force from the 
laser beams [1-3]; for modulated beams it periodically 
depends on time, G(i -I- tf) — G(i), with Tp = 2ti /ujp 
being the modulation period. The force {{t) in Eq. (1) 
is the zero-mean noise from spontaneous light emis- 
sion, which is often modeled by white Gaussian noise, 
(i')> = "S^makBTd^^^Sit - t'), where T is the 
effective friction coefficient, T is the characteristic gas 
temperature, and k,k' = 1,2,3 enumerate the atomic 
coordinates. 

We are interested in the parameter range where, in 

the absence of noise, Eq. (1) has two stable periodic so- 
lutions ri_2(i), with r„(i -f 2rj?) = r„(t) (n = 1,2) and 
r2{t) = ri{t -\- Tp)- These solutions give the positions 
of the centers of the vibrating atomic clouds in Fig. la 
of the main paper. By linearizing Eq. (1) about r„(t) 
one finds the typical frequencies ujci and relaxation rates 
(~ r) of the intracloud motion. Noise leads to fluctua- 
tions about r„ (t) with variance oc T that determines the 
size of the clouds in Fig. 1. 

Switching between the clouds results from large rare 
fluctuations. It is central for the analysis that the switch- 
ing rates are much smaller than ujf and than a;ci,r. 
Therefore of interest are period-averaged rates Wnm of 



n ^ m switching, n, m = 1, 2. Quite generally [4, 5] 

Wnm = Cw exp{-Rn/kBT), (2) 

where the prefactor Cyy ^ max(wci-r)- The most im- 
portant in Wnm is the exponential factor, which has the 
activation form. The activation energy ii„ of switching 
from the periodic state n is given by the solution of a 
variational problem i?„ = min??,, where the functional 
TZ is [5], 

TZ = (SFwa)"^ j dti{tf + j dtX{t) [mr - G - f{t)] (3) 

Here, the term quadratic in f determines the probability 
density of a realization of the random force that leads to 
switching [6]; X{t) is the Lagrange multiplier that relates 
the dynamics of the atom and the force to each other. 
For switching from nth periodic state, the trajectory r{t) 
should go from r„(t) for t — > — oo to the unstable periodic 
saddle- type state r(,(f) on the boundary of the basin of 
attraction to r„(t) for f — >■ oo; X{t), f (t) — >■ for t — >■ ±oo. 

Since in the single-atom case the period-two vibra- 
tional states difi'er only by phase, the activation ener- 
gies R\^2 are equal, R\ = R2 = R'^^\ The variational 
problem is significantly simplified if the vibrations r„(f) 
are almost sinusoidal, which corresponds to small damp- 
ing and comparatively weak nonlinearity, see below. For 
a one-dimensional weakly nonlinear oscillator was 
found and measured earlier [7-9]. 

B. Switching rate modification by the atom-atom 
interaction 

The interaction between atoms leads to an extra force 
Giiit(i") in Eq. (1). This force is weak in the sense that 
it weakly affects the intracloud atomic dynamics. In the 
experiment, the radius of the clouds changed by ^ 10% 
where the number of atoms changed by a factor of 2. 
The change of the vibration amplitude was also small, 
see Fig. 1 of the main paper. 

The force Gint depends on the positions of all atoms. 
In the analysis of the effect of the interaction on the in- 
tracloud dynamics it can be separated into the mean- 
field and fluctuating parts. The mean-field part is due 
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to the long-range interaction and is described in terms of 
the continuous atomic density [1]. The fluctuating part 
is dominated by short-range elastic collisions. For com- 
paratively low atomic densities and high temperatures 
used in the experiment the random part of Gjnt is small 
compared to the thermal noise and will be disregarded. 
The inter-cloud interaction is determined by long-range 
forces and can be described in the mean-field approxima- 
tion: the force Gi„t(r) from another cloud changes only 
slightly when atoms in this cloud shift by an interatomic 
distance. 

From the above arguments, the interatomic force on 
an atom in cloud n, Gint(r; n), can be written as 

Gint(r; n) = iV^Gintlr; n|m), (4) 

where Nm is the number of atoms in cloud m. The force 

Gint (r;n|m) is determined by the density distribution in 
the cloud m, and the terms with n = m and n ^ m 
describe the force from the atoms in the same and in the 
different cloud, respectively. 

Even though Gi„t is small compared to the single- 
atom force G, it may strongly affect the switching rate, 
leading to the spontaneous symmetry breaking. This 
is because Gint changes the switching activation energy 
i?„. Extending to the many-atom system the approach 
of Ref. 5, to first order in the interaction we obtain 
Rn = R^°^ + with 

m 

/oo 
dal°P')(t)Gi„t(r(,°P*)(0;n|m). (5) 
-oo 

Here, r^'^^\t) and \^n'^*'\t) describe the trajectory that 
minimizes functional 7^, Eq. (3), for switching from state 

n. Function rn'^^\t) gives the most probable path fol- 
lowed in switching; for a parametrically excited oscillator 
this path was recently observed in experiment [10]. Since 
the states 1 and 2 differ only in phase, we have an = a22 
and = 021- 

The following arguments were used in deriving Eq. (5) . 
First, because the interatomic interaction is weak, atoms 
switch one by one. not in groups; switchings of individual 
atoms are uncorrelated. Second, and most importantly, 
in the optimal fluctuation leading to switching of an atom 
the atomic clouds practically do not change their shape. 
This is also because the interatomic interaction is weak. 



and the force from a single atom in the cloud very weakly 
affects the motion of the switching atom, jormil ^ ksT. 
Only the effect from all atoms in the cloud is apprecia- 
ble, the product anmNm (-^1,2 ^ 1) can become larger 
than ksT. On the other hand, a significant change of 
the shapes of the clouds, even though it would change 
the switching activation energy, has an exceedingly small 
probability, which can not compensate the decrease of 
r'^iI- Indeed, if the cloud radius is ag ^ {ksT /rUaOjI^Y^'^ , 
the probability density to have the atoms in nth cloud 
displaced by Ar„ » ao is oc exp(— iV,iArJ;/2aQ). There- 
fore in calculating Gi„i (rn°''''(i); n|m) one should assume 
that the atoms in clouds 1 and 2 are at their stable equi- 
librium positions Vi{t) and r2(t). In other words, there is 
no fluctuational barrier preparation [11, 12] for intercloud 
switching, even though the switching barrier is strongly 
affected by the many-body interaction. This significantly 
simplifies the explicit calculation. 

Equation (5) shows that the effective switching activa- 
tion energy linearly depends on the numbers of atoms in 
the clouds, for weak interatomic coupling. Even though 

|i?n ^ I is small compared to the single-atom activation en- 
ergy R^^\ it can largely exceed fc^T, in which case the 
change of the switching rate will be large. From Eqs. (2) 
and (5) we obtain for the many-atom switching rate 

A^tot - A^i + A^2, (6) 

where = Cw exp(— /T) is the switching rate in 
the neglect of the interatomic interaction; a = (an — 
ai2)/2kBT and /3 = -(an + ai2)/2kBT. Equation (6) 
coincides with Eq. (1) of the paper and provides a general 
form of the coefficients a,j3. 



II. THE SYMMETRY BREAKING 
TRANSITION 

A. The master equation 

Because atoms move fast within the clouds compared 
to inter-cloud switching, on time scale ~ the state 
of the system is fully characterized by the mimbcrs of 
atoms in each cloud Ni,N2. For a given A^tot) the system 
is described by probability Pi{Ni) to have A''i atoms in 
cloud 1 at time t. This probability can be found from 
the master equation 



A(^i) = -[^(iVi) + iy{N,)]Pi{Ni) + iy{Ni - l)Pi(iVi - 1) + ^i{N, + l)Pi{N, + 1), 

f,(Ni) = NiWi2{Ni;N,ot), i^{Ni) = (N.^t - A^i)W^i2(A^tot - A^i; A^tot)- (7) 

I 

Here we have taken into account that, as a result of a creases by 1 and in cloud m increases by 1, and that the 
n ^ m transition, the number of atoms in cloud n de- 
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total probability of such a transition per unit time for 7V„ 
atoms in cloud n is NnW„m- 

Prior to making a transition the atom will move ran- 
domly within the cloud for a long time and will com- 
pletely "forget" its initial state. It is this randomization 
that makes the many-particle system zero-dimensional 
and leads to the mean- field approximation being exact. 

The stationary solution of Eq. (7) is 

Pf (TVi) = Z-^ (^^°*^ exp [-2aN, {N,,, -N,)], (8) 

where Z = Z{Ntot) is the normalization constant given 
by condition J^Ni Pi{Ni) = 1. 

For large A^tot, the function Pl^{Ni) has one or two 
sharp peaks. The location of the peak(s) depends on the 
total number of atoms. We will study the critical region 
where A^tot is close to 1/a, in which case near the peak(s) 
|A^i — iV2| <C iVtot- It is convenient then to introduce a 
quasi-continuous variable x = {N2 — Ni)/Niot = 1 — 
2Ni/Ntof For |a;| < 1 from Eq. (8) we find 

Pf (TVi) « Z-i exp [-Ntotix^ - 6fe2)/i2] , 

e = aN,ot - 1, (9) 

where Z is the normalization constant. 

Equation (9) has the standard form of the mean-field 
probability distribution near a symmetry-breaking tran- 
sition [13]. The parameter 6 is the control parameter, 
which plays the role of the deviation from the critical 
temperature in systems in thermal equilibrium, whereas 
X plays the role of the ordc;r parameter. Parameter 6 lin- 
early depends on the total number of atoms; it also lin- 
early depends on the reciprocal temperature. The critical 
number of atoms Nc is determined by condition ^ = 0, 

Nc = l/a(xT. (10) 



B. Critical exponents 

The mean reduced difference of the cloud populations 
T] = (x) is determined by the position xq of the maximum 
of the distribution Pf{Ni). For 6 < the distribution 
has one maximum at Xq = t] = 0, whereas at > it 
has two symmetric maxima. The system occupies one of 
them, which corresponds to the spontaneous symmetry 
breaking. Close to the critical point 

r; = Ar-i(iV2-7Vi) = ±(361)1/2 for < 6I < 1, (11) 

which is the familiar mean-field scaling with the critical 
exponent 1/2. 

Close, but not too close to the critical point, 1 ^ 

1 /2 

\0\ > Nc ' , from Eq. (9) we obtain for the mean 
square fluctuations of the populations difference cr^ = 
N-J{{N2-N,f)-rj' 

= (7V,|6l|)-i for 6* < 0; 
a^ = {2Ne\0\y^ iore>0 (12) 



(1^1 <C 1). This shows the familiar scaling of the 
variance of the order parameter on the both sides of 
the critical point. At the critical point, 9 = 0, we have 
cr2 = 27rV6N~^^^/Til/4f » 1.2iV^r^^^ Note that (x^) 
remains finite in a finite-size system at the critical point, 

— 1/2 

but its dependence on A'tot ~ -^c is given by factor N^^^ ' 
instead of N^l far from criticality. 

C. Response to the symmetry-breEiking field 

The symmetry of the period-two vibrational states can 
be lifted if the system is driven by an extra additive force 
h{t) = hcos(wFV2 + 4'h)- Such force is an analog of a 
static force in a thermal equilibrium system, it is "static" 
in the frame oscillating at frequency u)f/2. In the experi- 
ment the force was produced by modulating the counter- 
propagating laser beams at frequency ujf/2 with different 
amplitude; the phase <j)h is counted off from the phase of 
the strong modulation with period Tp. Response to a 
period-2rj' force becomes strongly nonlinear at the criti- 
cal point 9 = even for small amplitude |h|, because the 
force modifies the switching activation energies i?i,2 [14]. 
The latter can be found from Eq. (3) with G G + h(t). 
Extending the analysis of Refs. 5, 14, we obtain to first 
order in the inter-atomic interaction and the symmetry- 
lifting field 

P„ = iiW + iiW+pW, 

= - j da(,°P*' (t)h(t) = h COS (j)n. (13) 

The terms R^^\ are proportional to the field amplitude, 
h (X |h|. The phases (/)i^2 are the same as for a single 
parametrically modulated oscillator additionally driven 
by a force at frequency wf/2 [14]. They are linear in (t>h, 
with 02 = + TT. This can be immediately seen from 
Eq. (13), since by symmetry A^°P*^(t) = A^°P*'(t tf), 
whereas h(i + tf) = — h(i)- As a consequence, R^i"^ = 
—R^\ The difference in the activation energies of 1 ^ 2 
and 2 — >■ 1 transitions lifts the degeneracy of the states 
1 and 2, making their populations different even in the 
absence of the inter-atomic coupling. 

Even though \R^n \ < for small field amplitude 

|h|, the ratio /UbT should not be small, and there- 
fore the field-induced change of the switching rates may 
be substantial. At the same time, the change of the am- 
plitudes of cloud vibrations remains small for small h and 
will not be discussed. Equation (13) for R^rf* resembles 
the change of the free energy of an Ising spin by a mag- 
netic field oc h tilted by an angle 0i with respect to the 
quantization axis. 

Formally, the terms R^^^ lead, respectively, to fac- 
tors exp[—(/i/fcBT) cos0i] and exp[(/i/fcB7') cos0i] in the 
rate coefficients n{Ni) and i^iNi) in master equation (7). 
As a result the stationary distribution Pf'(7Vi), Eq. (9), 
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acquires an extra factor 

Pf (ATi ) ^ Pf (TVi ) exp [- (Ntothx/kBT) cos .^i] (14) 

[the; normalization constant in Pj** is changed appropri- 
ately; we remind that x = {N2 — Ni)/Niot-]- 

The equation for the maximum of the distribution 

xq = T] near criticality now becomes 9xo — -^Xq — 
(h/kBT) cos4>i = 0. For small \h/kBT\ < 1 and not too 

— 1/2 

close to the phase-transition point, 1 \9\'^ Nc , we 
have 

rj K 6-'^ {h/kBT) cos (t>i {e<0), (15) 
T) « ±(36»)V2 _ (26)-'^ (h/kBT) cos?!>i {6 > 0). 

The generalized resonant susceptibility S(0) = —drj/dh 
scales as \9\~^ on the both sides of the phase transition 
point and diverges for ^ — >■ 0. We note that, even though 
wc consider a system far from equilibrium and the sus- 
ceptibility is calculated with respect to a time- dependent 
field, for A/tot = const it is simply related to the variance 
of the order parameter cr^ as given by Eq. (12). 

At the phase transition point, 9 = Q, there is no linear 
static susceptibility, and rj = xq <x h}/^, as seen in the 
experiment. 

D. Frequency dispersion of resonant response 

The above analysis can be extended also to the case of 
an extra additive force with frequency (wf/2) -|- O that 
differs from, but remains close to wf/Z, i.e., |n| <C F ^ 
ujp. To find the distribution Pi near the maximum, one 
can transform Eq. (7) into a Fokker-Planck equation by 
expanding Pi(7Vi ± 1) - Pi{Ni) « ±dN,Pi + (1/2)9^^ Pi 
and using similar expansions for /z, v [15]. In the absence 
of extra modulation, for small x = {N2 — Ni)/N^ot to 
leading order in x the equation reads 

W-^Pi = LPu W = exp(^Artot), 

LPi = -2d^ [{9x - lx^)Pi] + 2N,-ldlPi. (16) 

The stationary solution of Eq. (16) Pf* has the form (9). 

An extra additive force hcos [(^0;;- -t- 17) with small 
fl adiabatically modulates the rates Wnm of the inter- 
cloud transitions, i.e., one can think of the rates Wnm as 
parametrically dependent on time. The major contribu- 
tion to the rate modulation comes from the modulation 
of the activation energies, which is described by Eq. (13) 
with 

r[''^ = R^^^ (t) = -Rf^ (t) = h cos(Oi) (17) 

(we have disregarded a time-independent term in the ar- 
gument of the cosine; for quantities that do not oscillate 
at frequency ujp/'^, like cloud populations, this term can 
be eliminated by changing time origin). 

For weak modulation one can linearize the transition 
rates in h and consider the terms oc /i as a perturba- 
tion. Then the distribution can be sought in the form 



Pi = Pi* -I- (5 Pi, where the field-induced term SPi oc /i is 
oscillating at frequency f2. It is given by equation 

W-^6Pi - LdPi 

= 2^ cos(Qt)Artot (^9x - ^x^^ Pf . (18) 

Equation (18) allows one to find the periodically oscil- 
lating term in the order parameter Sr]{t) = {x) — xo {xq is 
the position of the maximum of the stationary distribu- 
tion) and the generalized susceptibility !B(f2), which can 
be defined by expression 

Sriit) = [E{fl)hexp{-iat) +C.C.] . 

This can be done by multiplying Eq. (18) by x and inte- 
grating over X. Expanding 9x — x^/3 in L about xq, for 
1 > |6'| > iVt^y^ one finds 

E{n) ^2W/{2\9\W -in), 9<0, 

E{n) = 2W/i49W -in), 9>0. (19) 

[in the critical region, in Eq. (16) for W one should re- 
place A^tot with Nc]. For O / the suscoptibihty Eq. (19) 
remains finite even at the phase transition, 9 = 0. Here 
E{n) = 2iW/n. The susceptibility diverges with fre- 
quency, with critical exponent equal to —1. 



III. FLUCTUATIONS OF THE TOTAL TRAP 
POPULATION 

In the above analysis we disregarded inelastic collisions 

between trapped atoms and collisions with atoms out- 
side the trap. These processes lead to fluctuations of 
the total number of atoms in the trap. In the exper- 
imental conditions they were slow, with rates an order 
of magnitude smaller than the rate of inter-cloud transi- 
tions Wnm- However, near the critical point inter-cloud 
population fluctuations are slowed down, and then fluc- 
tuations of the total number of trapped atoms A^tot may 
become important. 

In the presence of fluctuations of A'^tot the state of the 
system should be characterized by the probability P(N) 
[N = {Ni,N2)] to have A'l atoms in cloud 1 and N2 
atoms in cloud 2, with A'tot = A^i -I- A''2. For compara- 
tively small A^tot and low temperatures used in the experi- 
ment, evaporation from the trap due to inelastic collisions 
between trapped atoms was rare. The major mechanism 
of the change of A^tot was collisions with atoms oiitside 
the trap. In a simple model of such collisions, the rate of 
loss of atoms from an nth cloud Wout-^n is proportional 
to the number of atoms in the cloud Nn, whereas the 
rate of capture of outside atoms by a cloud Wout-^o/2 
is independent of A^„ and depends only on the density 
of outside atoms (which determines A'o) and the cross- 
section of their capture. We disregard fluctuations in 
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the density of outside atoms and choose parameter Nq in will be seen below. 

such a way that it provides the typical scale of iVtot, as The master equation for F(N) has the form 



P(N) = WLaP + WontLontP, LeiP = -V /x'(7V„;7Vtot)P(N) 

^ — 

+ 1; ^tot)P(iVi + 1, 7V2 - 1) + fi'iN2 + 1; Ntot)P{Ni - 1, iVa + 1), (20) 
LoutP = -(iVi + iV2 + iVo)P(N) + V (iV„ + l)P(N + 5„) + (iVo/2)V P(N-5„), 



where 

n'{Nn;Ntot) = iV„exp {aNtot - ^aN^) = W^-V(^«) 

and ^1 = (1, 0) and ^2 = (0, 1) [the inter-cloud switching 
rate W is defined in Eqs. (6), (16)]. The terms L^x and 
-^out describe inter-cloud transitions and exchange with 
atoms outside the trap, respectively. We assume that the 
typical rate of this exchange Wout W. 

We now consider the critical region where A^tot is close 
to Nc- We will be interested in the range of iVi, A^2 where 
P{Ni,N2) is close to maximum, \N-i — <^ -^tot, and 
introduce quasi-continuous variables 



X = 



N2-N1 

No ' 
aNo - 1. 



u= vj 1, (21) 

iVo 



As we will see, variables x and 9 coincide with the vari- 
ables used earlier if fluctuations of the total number of 
trapped atoms can be disregarded; variable u gives the 
deviation of iVtot from Nq. 

For \x\ <^ 1, 1^1 <^ 1, and \u\ ^ 1 to leading order in 
x,6,u, and Nq^ the operators L^i and Lo^t become 



LoiP 



LniitP 



-2d, 



x{9 + u)- Ix^ 



+2N^'dlP, 



(22) 



d^xP) + d^iuP) + N^' [{dl + dl)p] . 



Equations (20) - (22) allow one to study the stationary 
distribution of the modulated system in the presence of 
fluctuations of the total population. We note first that 
exchange of atoms with the surrounding leads to effec- 
tive exchange of atoms between the clouds. In turn, this 
leads to rcnormalization of inter-cloud transition rates 
and, effectively, of the control parameter, 



6^0 



2W 



a 



1 + 



2W 



(23) 



The change of 9 is negative, which should be expected, 
since exchange with the surrounding should stabilize the 

symmetric phase, and therefore a larger number of atoms 
is required for the symmetry breaking transition than in 



the case where there is no such exchange. Nc is now the 

critical value of the parameter Nq, the mean number of 
trapped atoms is {Ntot) = Nq{1 + {u)). The coefficient of 
diffusion over x is also renormalized, 2Nq^ — >■ 2Nq^{1 + 

Wont/2W). 

In the stationary regime fluctuations of the total num- 
ber of atoms in the trap are small, (u^) ~ ^/Nq. How- 
ever, the distribution over the scaled difi^erence of the 
cloud populations x differs from the standard Landau- 
type distribution (9). This is in spite the fact that the 
system is described by the mean-field theory, no spatial 
correlations are involved. The dynamics of x and u are 
very different. There is no critical slowing down for u, 



and in the critical region \9\ < N 
distribution over x is ^ N, 



'1/2 



the width of the 



A. The variance of the order parameter 



The major effect of the coupling between fluctuations 
of X and u is the change of the variance = (x^) — (x)^. 
The latter occurs in the broken-symmetry state where 
9^Nc'^''. It can be calculated in the stationary regime 
by decoupling the chain of equations for the moments of 
the stationary distribution {x'^u'') with integer k,l. For 



-1/2 



<6'< 1 



1 3 

29 49 + (Wont /W)' 



(24) 



Fluctuations of u and x (i.e., of the total population 
of the clouds and the population difference) are corre- 
lated in the broken-symmetry phase, with {u{x — xq)) « 
{2xo/N,)[49 + {Wont/W)]-\ where xq = ±{39)'^^. 
In contrast, in the symmetric phase for 1 ;:g> — ^ 



N- 

Nr 



-1/2 



fluctuations of u and 



\9\- 



1 



are uncorrelated, 

.\9\-\ 



with 



(Wout/2VK) 

It is seen from Eq. (24) that, for 9 > Wont/4W, 
the variance of the order parameter fluctuations in the 

broken-symmetry state scales with the control parame- 
ter in the same way as without fluctuations of the total 
population, Nca^ w 5/4:9. However, the factor in front 
of 1/9 is now larger than in the symmetric phase. This 
corresponds to what was seen in the experiment. As 9 
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decreases below Wout/4W^, there occurs a crossover to 
« l/2Nc9, i.e., scales with 9 with the same scaling 
exponent, but with smaller amplitude than in the sym- 
metric phase. 

We note that there are other fluctuations that can con- 
tribute to the variance cr^. For example, slow fluctuations 
of the trap parameters may lead to small and slow fluc- 
tuations of the total population, which can be thought 
of as fliictuations of the control parameter 9. They give 
an extra contribution 3{{59)^) /M.9) to in Eq. (12) for 
the symmetry-broken phase, where {69'^) is the variance 
of 9. Interestingly, this contribution also scales as 9'^ 
away from the critical point. It further increases the dif- 
ference between the variance of the order parameter in 
the broken-symmetry and symmetric phase. 



IV. ONE-DIMENSIONAL MODEL 



A. The rotating wave approximation 

Calculation of the switching rates and response coef- 
ficients is simplified if the atomic motion in MOT is as- 
sumed one-dimensional, along the MOT axis z. This cor- 
responds to averaging over the motion transverse to the 
axis, which is largely small-amplitude fluctuations about 
the vibrating centers of the clouds. In the experiment 
the radiation that confined the atoms transverse to the 
axis had large intensity, so that the correlation time of 
these fluctuations was short. In the single- atom approxi- 
mation the atomic motion can be modeled by the Duffing 
oscillator, 

z + 2Tz + Wo(l + ep cosujptjz + -fz^ = f{t)/ma, (25) 

Here, coq and F are the MOT eigenfrequency and viscous 
friction coefficient, respectively, 7 is the nonlinearity pa- 
rameter, and ep is determined by the amplitude of the 
laser beam modulation [16]. These parameters are given 
in the main paper. The terms nonlinear in the veloc- 
ity, which are disregarded in Eq. (25), arc comparatively 
small for the experimental conditions. The function f{t) 
is white thermal Gaussian noise, cf. Eq. (1). 

For the reduced modulation strength epUJo/Ar > 1, 
resonant modulation (jwi? — 2a;o| ^ wq) can excite period- 
two atomic vibrations [17]. The vibrations are nearly si- 
nusoidal in the absence of noise, if the modulation is not 
too strong. The atomic dynamics can be conveniently de- 
scribed by switching to the rotating frame using a stan- 
dard transformation 

z = Crwa [92 cos(w_Ft/2) - qi sin(w_Ft/2)] , (26) 
i = -Crwa(wf/2) [q2 sm{uiFt/2) + qi cos{uiFt/2)] 

with Crwa = (2wQei?/37)^/^ (the modulation sign is cho- 
sen in such a way that ep/j > 0). In the rotating wave 



approximation (RWA), the equations for slow variables 
q = (gi, 92) in slow time r are 



rfq 



= K(q) + f(T) 



(27) 



Here, f(T) is white Gaussian noise with two asymp- 
totically independent components, {f k{t) f n' {t')) = 
2DrkBT6^^'5{T - r'), with k,k' = 1,2 and Dr = 
Q^V /rriaOjQe^. We use vector notations here formally, 
q is not a vector in real space, its components arc com- 
binations of atomic displacement and momentum along 
the MOT axis. 

Functions K(q) are cubic polynomials and are given 
explicitly in Ref. 14 [they were denoted by K(0)(q)]. The 
zeros of K occur at the stationary states in the rotating 
frame and, in turn, at the period-two vibrational states in 
the laboratory frame. We will be interested in the param- 
eter range where Eq. (27) has two attractors q^ = — q^. 
Functions K are also equal to zero for q = 0, which corre- 
sponds to the unstable state of zero-amplitude period-two 
vibrations in the laboratory frame. 

For low temperatures, the motion described by 
Erj. (27) is primarily small-amplitude fluctuations about 
states q^2- Switching between the states requires a large 
rare fluctuation. The single-atom switching rate has ac- 
tivation form T4^(") = Cexp(-i?(")/fcsT), Eq. (2), where 
as in Eq. (3), = min 7^(°) with 



7e(°) = {4Dr)-^ J rfrf2(r) 
dq 



I 



+ / drA(r) 



dr 



K-f(r) 



(28) 



The major distinction from the general formulation (3) 
is that function K does not explicitly depend on time. 
This made it possible to solve the variational problem 
(28) and to find both the optimal trajectories followed in 

switching from nth stable state qn°^*^ (r) [7] and functions 
A„(r), which are equal to the logarithmic susceptibilities 
of Ref. 14 multiplied by -l/Dr- 



B. Interatomic interaction 

The major inter-atomic interaction that affects the 

switching rate is the long-range interaction, since in 
switching atoms move far away from the clouds. It comes 
from the shadow effect, which is due to "shielding" of 
atoms from the laser light by other atoms [1, 2, 18-20]. 
In the one-dimensional picture, the force on ith atom 
with coordinate from jth atom with coordinate z^ can 
be approximately modeled as = —fsh J2j sgn(2;* — z^). 
This force is weak, much smaller than the Doppler force 
that confines atoms to the trap. Multiple scattering of 
light is disregarded in the above expression. 

To estimate fsh for a switching atom we have to take 
into account that the atomic clouds are 3-dimcnsional. 
We consider a simple model in which a laser beam with 
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intensity / propagating along z-axis passes on its way 
through an atomic cloud with density distribution p(r). 
The resulting intensity change as a function of transverse 
coordinates x,y is AI{x,y) = Igl J dzp{r), where inte- 
gration is done over the length of the cloud and is 
the absorption cross-section. This cross-section depends 
in the standard way on the light intensity and the fre- 
quency detuning. Generally, because of the magnetic- 
field induced frequency shift and the Doppler shift (t^ 
oscillates in time. The light intensity / also oscillates in 
time. However, for comparatively small modulation am- 
plitudes, we can disregard these oscillations. Then we get 
for typical experimental conditions ~ 5.6 x 10~^^ m^. 

The atomic density distribution p{r) can be assumed 
Gaussian with the same width Wt ^ 1 mm in all direc- 
tions, which was achieved in the experiment by tuning 
the transverse beam intensities. According to the op- 
timal path picture, during switching atoms most likely 
move along the MOT axis, and for such atoms the light 
intensity change is determined by p{r) on the axis. The 
extra force on the switching atom as it moves between 
the clouds is then directed toward the more populated 
cloud and is equal to /sh|-^i — -^2!, with 

U = hkTpaLs/4TT{s + l)wl (29) 

Here, k is the photon wave number, Tp is the recipro- 
cal lifetime of the excited state, and s = {I/Is)[l + 
{26/Tp)'^]~^ is the resonant absorption strength {Is is the 
saturation intensity and 6 is the detuning of the radiation 
freqiicncy from the atomic transition frequency). Equa- 
tion (29) gives » 2.5 X 10"^^ N. Note that, since 
the clouds are oscillating, the force oc fsh is oscillating 
as well, its time dependence is determined by the sgn- 
function. 



C. The shadow effect in the rotating frame 

The effect of the shadow force can be conveniently an- 
alyzed in the rotating frame by changing from (2;*, i') to 
slowly varying two-component vectors q', Eq. (26). In 
the RWA the equation of motion for q* becomes 

= K'(q') + ed^iH,^ + r(T) (30) 

{i — 1, . . . , iVtot)- Here, f*(r) is the random force on 
ith atom (random forces on different atoms are uncorre- 
lated), e is the permutation tensor, and iJgh is the inter- 
action Hamiltonian in the slow variables, 

^^sh=^V'..V;h(q'-q^), 

l^sh(q) = |q|, (31) 



where the prime indicates that i ^ j. 

The general expression for the interaction-induced cor- 
rection to the activation energy ii„, Eq. (5), in the RWA 
has the form 

/oo 
rfrA„(T)eSq^y«h(qi°P*n^) - q^), (32) 
-CX) 

where A„ (r) and qi°^''' (r) are the solutions of the single- 
atom problem of minimizing the functional (28) for 
switching from the state n. 

Equations (28) and (32) were used to find the coef- 
ficients a„„i and the critical value of the total number 
of atoms Nc = 2kBT / {an — ai2) where the symmetry 
breaking transition occurs. The obtained Nc was within 
a factor of 2 from the value of Nc observed in the exper- 
iment. This is reasonable given the uncertainty in the 
temperature, the renormalization of Nc by the finite life- 
time of atoms in the MOT, cf. Eq. (23), and the fact 
that the light intensity modulation was not weak and, 
respectively, the atomic vibrations noticeably deviated 
from sinusoidal. 

An important effect that allows one to independently 
compare the many-atomic theory with experiment is the 
change of the vibration amplitudes due to the interac- 
tion. This change is determined by the change of the 
equilibrium positions q,^ in the rotating frame. For weak 
inter-atomic interaction it can be found by linearizing 
equations of motion (30). The average force on an atom 
from other atoms in the same cloud is equal to zero. The 
shift of the equilibrium position is due only to the force 
from the atoms in the other cloud. Disregarding small 
fluctuations about the equilibrium positions, for the shift 
of nth attractor in the rotating frame 8(\^ we obtain 

(<5q;tSq.)K {ai) + N^_^ed^AV,^{cit - q^ J = 0. (33) 

Here, Vsh is given by Eq. (31) for Vsh. in which /sh is 
replaced with /sh. The quantity /sh characterizes the 
force from the shadow effect which is averaged over the 
cross-section of the cloud transverse to the MOT axis. 
Such averaging occurs as atoms arc moving within the 
cjoud. For Gaussian clouds of the same cross-section, 

/sh = /sh/2. 

The shift of the nth attractor (n = 1,2) given by 
Eq. (33) and the corresponding change of the vibration 
amplitude of the nth cloud is proportional to the num- 
ber of atoms in the other cloud, A^3_„. Therefore the 
vibration amplitudes of the clouds become different in 
the broken-symmetry state. This prediction and the lin- 
ear dependence of the amplitude change on the number 
of atoms in the other cloud are in full agreement with 
the experiment; the numerical estimate for the simplified 
model of sinusoidal ID vibrations is within a factor of 2 
from the measured value. 
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